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We investigated the photocurrent in poly(3-hexylthiophene-2,5-diyl) (P3HT):[6,6]-phenyl-C6i butyric acid 
methyl ester (PCBM) solar cells by applying a pulsed measurement technique. For annealed samples, a point 
of optimal symmetry (POS) with a corresponding voltage Vpos °f 0.52-0.64 V could be determined. Based on 
macroscopic simulations and results from capacitance-voltage measurements, we identify this voltage with flat 
band conditions in the bulk of the cell, but not the built-in voltage as proposed by [Ooi et al., J. Mater. Chem. 
18 (2008) 1644]. We calculated the field dependent polaron pair dissociation after Onsager-Braun and the 
voltage dependent extraction of charge carriers after Sokel and Hughes with respect to this point of symmetry. 
Our analysis allows to explain the experimental photocurrent in both forward and reverse directions. Also, we 
observed a voltage-independent offset of the photocurrent. As this offset is crucial for the device performance, 
we investigated its dependence on cathode material and thermal treatment. From our considerations we gain 
new insight into the photocurrent 's voltage dependence and the limitations of device efficiency. 

PACS numbers: 71.23.An, 72.20.Jv, 72.80.Le, 73.50.Pz, 73.63.Bd 
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I. INTRODUCTION 

Organic solar cells have improved greatly in the last years, 
reaching 5-6 % power conversion efficiency (PCE) today. 1 
However, the Shockley diode equation cannot explain the 
voltage dependent photocurrent in organic solar cells based 
on a physical model™ The detailed process leading from pho- 
toinduced polaron pairs to extracted charge carriers and ex- 
ternal photocurrent needs better understanding to push devel- 
opment of these devices systematically. 3 The first step in this 
process is singlet exciton formation upon absorption of a pho- 
ton, usually on the polymer. Then, the exciton has to diffuse 
to a polymer-fullerene interface within its lifetime, where, due 
to the high electron affinity of the fullerene, a polaron pair is 
created via ultrafast charge-transfer. 4 These polaron pairs are 
still Coulomb-bound due to a low dielectric constant £,- of the 
organic material system (typically 3^1), and have to be disso- 
ciated to free polarons. Provided the free polarons arrive at 
their respective electrodes, the last step is charge extraction. 

In this paper, we investigate the symmetry and voltage de- 
pendence of the photocurrent, and compare this to a model 
that takes field dependent polaron pair dissociation and volt- 
age dependent charge extraction into account. The photocur- 
rent Jph, defined as the difference of illuminated and dark cur- 
rent, /ught — ./dark* can experimentally be accessed using pulsed 
illumination (inset in Figure [IJ. This is necessary to avoid a 
major overestimation due to device heatingP 

We find the photocurrent of annealed P3HT:PCBM solar 
cells to have a point of optimal symmetry (POS). The pho- 
tocurrent shows a point symmetry with respect to POS (see 



Figure [TJ, and the corresponding voltage Vpos is in the range 
of 0.52-0.64 V. While these values agree with results of Ooi 
et alP\ our interpretation differs. We identify this voltage as 
the quasi flat band case-with flat bands in the bulk of the so- 
lar cell-which is well below the built-in voltage. Accordingly, 
our considerations lead to new insights into the origin of the 
photocurrent. 

As shown in Figure [T] the photocurrent is composed of two 
contributions. The first is voltage-dependent and symmetric 
with respect to POS. The second contribution is a voltage- 
independent offset. This constant offset is usually negative, 
thereby increasing the short-circuit current 7sc an d PCE. As 
the magnitude of this offset is crucial for device performance, 
we investigated its dependence on thermal treatment, cathode 
material and active layer thickness in detail. 

In order to describe the voltage dependence of 7ph, we use 
a combination of Braun-Onsagei 6 7 theory and charge extrac- 
tion as calculated by Sokel and HughesP This combination 
was proposed by Mihailetchi et al. 9 and could for the first 
time explain the experimental photocurrent in reverse direc- 
tion. However, they calculated effective voltage and field 
with respect to the physically ill-defined voltage Vb, at which 
/ph = 0. As we will show, this is correct only in the special 
case where Vo coincides with Vpos- In general, contact effects 
result in a lower value of Vpos- Consequently, we consider 
voltage and electric field relative to Vpos instead of Vb. Cal- 
culating polaron pair dissociation and charge extraction with 
respect to the point of symmetry now allows to describe the 
experimental photocurrent for solar cells with different elec- 
trode material. Using a pulsed measurement technique also 
makes the photocurrent in forward direction accessible, and 
shows that the model can describe the data in both directions. 
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FIG. 1. (Color online) Photocurrent density Jp^ vs. applied volt- 
age of a typical P3HT:PCBM solar cells with Ca/Al cathode under 1 
sun. The point of optimal symmetry (POS) is indicated by the blue 
circle. The voltage-dependent, symmetric shape of Jph is shifted by 
a constant offset to more negative values. Note that the density of 
recorded data points is higher between V and 1 V applied voltage. 
The inset shows the current flowing under pulsed illumination, with 
an increase in the applied voltage every millisecond. 



II. EXPERIMENTAL 

All samples investigated were spin-coated from solution 
of a 1:1 weight ratio of P3HT and PCBM in chlorobenzene, 
onto indium tin oxide/glass substrates coated with poly(3,4- 
ethylenedioxythiophene):poly(styrenesulfonate). After an an- 
nealing step at 140 °C, metal contacts were applied by ther- 
mal evaporation; either Ca (3 nm) followed by Al (1 10 nm) or 
only Ag (120 nm). All fabrication and characterization took 
place in a nitrogen atmosphere glovebox. The cells with Ca/Al 
cathode had fill-factors (FF) of about 60 % and PCE of 2-3 %, 
as determined with a Xe-lamp, adjusted to simulate standard 
testing conditions (STC).^ 

The photocurrent measurements were performed based on 
the setup proposed by Ooi et alR using a white 10 W LED 
for pulsed illumination. The pulse duration was set to 0.5 ms 
at a duty cycle of 50 %, giving the cell ample time to reach 
steady-state after switching the light on/off. The illumination 
level was equivalent to one sun, based on comparison of the 
short-cuircuit current (Jsc) reached with the LED to the 7sc 
reached under simulated STC. The active layer thickness was 
determined with a Veeco Dektak 150 profilometer. 

For the capacitance-voltage (CV) measurements we ap- 
plied a 5 kHz AC voltage with an amplitude of 40 mV to 
our cells and performed a DC bias sweep. It was carried out 
under nitrogen atmosphere and without any incident light to 
the device to prevent distortions due to charge carrier gen- 
eration. Data acquisition was carried out with an Agilent 
E4980A LCR-meter in parallel RC circuit mode. The rea- 



sonability of this working mode has been proven by preceding 
impedance measurements that showed a clear semicircle (with 
a negligible offset impedance of about 10 H, corresponding 
to the series resistance R s ) in Cole-Cole representation for 
our cells. Since we have Schottky-like contacts at the metal- 
semiconductor interfaces of our devices, the value of Vfcv was 
determined by the intercept of a linear fit of the Mott-Schottky 
plot. 

All presented simulations were done solving the one dimen- 
sional differential equation system of Poisson's equation and 
the continuity equations for electron and holes in a numerical 
iterative approachj^^l implicitly assuming a band-like trans- 
port. A 100 nm thick bulk-heterojunction solar cell is calcu- 
lated at a temperature of T = 300 K containing electrons and 
holes with equal mobilities of /j„ = p. p = 10 cm 2/VsPFor 
the active material, the effective medium approach is used, 
leading to effective electrical bandgap of Eg = 1.1 eV be- 
tween the donor HOMO (highest occupied molecular orbital) 
and the acceptor LUMO (lowest unoccupied molecular or- 
bital) level and an effective dielectric constant of e r = 3.4.E] 
The effective charge carrier densities of LUMO and HOMO — 
describing the spatial charge carrier densities of electrons 
n(x) and holes p(x) by Fermi -Dirac statistics — were set to 
N c = N y = 10 20 cm~ 3 . The interfaces towards the electrodes 
are assumed to possess injection barriers of 0.1 eV without 
any surface recombination.™ 

Illumination is taken into account by a homogenous charge 
carrier generation rate G over the whole device. This simpli- 
fication avoids band bending due to local differences of the 
charge carrier densities caused by optical interference inside 
the sample. We use Langevin theory to describe bimolecular 
charge carrier recombination of free polarons to the ground 
stated The combination of charge carrier generation and re- 
combination leads to the following net generation rate: 

U(x)=G — (np-rii) (1) 

fcOtr 

with the elementary charge q, the dielectric constant £o, the 
Boltzmann constant kg and the intrinsic charge carrier density 
rii = V" N C N V exp (qE G /2k B T). 

The result of this system are spatially resolved values for 
the electrical potential and electron and hole densities in form 
of quasi-Fermi potentials. In order to create a band structure, 
the potentials have to be transformed into an energy level. 
Since experimental values of device work functions are not 
available, the electron conducting contact (here aluminum) is 
set to a typical value of -4.25 eVH 6 which determines the ab- 
solute value of all other potentials. 

In Figure[2]the resulting band structures are presented in the 
dark and under illumination for two important bias voltages. 
First, a voltage resulting in flat bands — and therefore zero 
electric field E, being the spatial derivative of the potential — 
in the middle of the sample, and second the built-in voltage 
for comparison. 
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FIG. 2. (Color online) Energy structures of a simulated organic solar 
cell in the dark (left) and under illumination of 1 sun (right). The 
simulated HOMO and LUMO levels (solid lines) and quasi Fermi- 
levels for electrons and holes (dashed lines) of a 100 nm thick cell are 
shown, with the ITO contact on the left and the metal contact on the 
right side, a) At Vpos (QFB) the bands are flat in the bulk, while the 
electric field is finite in the contact regions, b) At an applied voltage 
equal to the built-in potential Vbi the vacuum levels at the position 
of the electrodes match, but the internal electric field is nonzero over 
the whole device. 



III. RESULTS AND DISCUSSION 
A. Origin of Vpos 

The photocurrent shows a high symmetry with respect to 
POS, with a corresponding voltage Vpos in the range of 0.52- 
0.64 V for annealed solar cells with Ca/Al cathode under 1 
sun. This agrees with findings of Ooi et al.Pwho reported a 
Vpos of 0.58-0.60 V for cells with Al cathode. 

To understand this symmetry in /ph, we look at the energy 
bands in the solar cell. Due to the band bending at the con- 
tacts, a case of flat bands in the whole device does not exist. 
However, our macroscopic simulation shows flat band condi- 
tions with zero electric field in the bulk of the cell at an applied 
voltage of 0.66 V (Figure [2^). We call this the quasi flat band 
(QFB) case, which is well below the built-in voltage Vbi at 
0.90 V (Figure^). 

The reason why the bias must be reduced from Vbi in order 
to achieve flat bands in the bulk is the band bending at the elec- 
trodes, which is a consequence of the boundary conditions for 
the quasi-Fermi levels at the interfaces. Ideally, the electrodes 
inject charge carriers of one type into the organic semiconduc- 
tor and extract the other type. This generates a strong diffusive 
current in the contact region which is, in thermal equilibrium, 
accompanied by a high electric field in the opposite direction, 
especially next to the electrodes. As shown in Figure [2] the 
resulting band bending is almost independent of illumination. 
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FIG. 3. (Color online) Typical result for a capacitance-voltage mea- 
surement of an annealed P3HT:PCBM solar cell with Ca/Al cath- 
ode in the dark. Normalized capacitance C/A (left axis) and Mott- 
Schottky plot (right axis). A linear fit to the Mott-Schottky plot re- 
sults in a value of 0.59 V for the quasi flat band voltage. 



Only the quasi-Fermi levels are affected, indicating a higher 
charge carrier density under illumination. 

This reduced voltage was published by Kemerink et al., but 
incorrectly identified as Vo. 17 At QFB, the electric field is zero 
in the bulk of the solar cell and finite in vicinity of the contacts, 
while at Vbi the field is finite at every position in the cell. 

Capacitance-voltage measurements — conducted in the 
dark — further support the interpretation of Vpos being the 
quasi flat band case. The capacitance C = SQ/hV — 
bQ/b(Ed) (where Q is the electric charge) depends on the 
differential variation of the electric field 5E in the bulk, which 
approaches zero in this case whereby the capacitance ap- 
proaches infinity. An extrapolation of C~ 2 to zero (linear fit 
to the Mott-Schottky plot, see Figure [3j yields Vcv> which 
therefore corresponds to the quasi flat band voltage. This is 
a common technique to determine the flat band voltage in 
inorganic semicondutorsp^ but has already been applied to 
organic 

devicesPEol We found V cv to be 0.5-0.6 V, which 
agrees with the observed values for Vpos- 

The quasi flat band case is distinct and can explain the sym- 
metry of the photocurrent. Energy bands with zero slope 
present a special, symmetric situation, as a relative bias in 
either direction will cause the photogenerated carriers in the 
bulk to flow equally strong, but in different directions. In con- 
clusion of this, we propose Vpos to correspond to the case of 
flat bands in the bulk of the active layer, with finite electric 
field in vicinity of the contacts, caused by band bending due 
to injection barriers. Hence, the electric field in the bulk will 
be proportional to an effective voltage \V — Vpos I- In contrast 
to Ooi et al.P we do not assume Vpos to equal the built-in 
voltage Vbi, but to be considerably smaller. 
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offsets for thinner and thicker devices (not shown). 

As an alternative explanation for this constant offset, we 
propose the band bending in vicinity of the contacts, which — 
independently of applied bias — has only one direction and re- 
sults in an electric field high enough for efficient polaron pair 
dissociation (see Figure [2]). The contact regions would then 
give a constant contribution to /ph, which depends on the de- 
gree of band bending — and therefore the height of the injec- 
tions barriers — while the bulk region contributes to 7ph as a 
function of applied voltage and causes the symmetry. 

The exact nature of this offset and its dependence on the 
cathode material has to be further investigated. A highly neg- 
ative offset is crucial for device performance, and should be 
optimized, e.g. by choice of cathode material. 



FIG. 4. (Color online) Photocurrent density vs. applied voltage, (a) 
Influence of thermal treatment for solar cells with Ca/AI cathode, 
(b) Cells with Ca/AI and Ag cathode, in each case for two different 
active layer thicknesses. The magnitude of the voltage-independent 
offset can be influenced by thermal treatment and choice of cathode 
material. 



B. Voltage-independent offset of 7ph 

As mentioned above, the photocurrent is composed of two 
contributions. The first is voltage-dependent and symmetric 
with respect to POS. The second contribution is a voltage- 
independent offset, which is crucial for device performance 
and can be influenced by processing parameters. Figure [4] 
shows the influence of cathode material and annealing con- 
ditions on the photocurrent. While the shape of 7ph remains 
similar, the offset is critically dependent on both cathode ma- 
terial and thermal treatment. With respect to the pristine sam- 
ple, the offset decreases after thermal treatment at 80 or 200 
°C, but increases upon treatment at 140 °C (Figure [4^). Us- 
ing Ag as cathode material yields a much smaller offset than 
Ca/AI, while the shape of the photocurrent curves is very simi- 
lar (Figure[4j5). Vpos is also smaller for cells with Ag cathode, 
with values of about 0.45 V, which is a direct consequence 
of the higher injection barriers. Active layer thickness d influ- 
ences the magnitude of the offset as well. The highest negative 
offset, up to 4. 1 mA cm" 1 , was observed for solar cells with 
Ca/AI cathode, active layer thickness of 120 nm and thermal 
annealing at 140 °C. These cells had the highest 7sc an d PCE. 

As this voltage-independent offset can be influenced by the 
choice of cathode material without changing the shape of the 
photocurrent significantly, we assume the contact regions to 
be responsible for the offset. The voltage dependent part of 
7ph then corresponds to the photocurrent generated in the bulk 
of the solar cell. 

Ooi et al. explained this offset with self-selective 
electrodes,^ resulting in a constant diffusion current 
7ph, offset = —(eDriph/d), where e is the elementary charge, D 
the diffusivity and n p h the concentration of photogenerated 
charge carriers at the selective electrode. In contrast to this 
predicted d^ 1 dependence, we observed a maximum (nega- 
tive) offset for a thickness d of about 120 nm, with smaller 



C. Voltage vs. field dependence 

The active layer thickness d has only little influence on 
the voltage dependence of the observed photocurrent (Fig- 
ure |5|a)). Of course the absolute value of /p n is higher for 
thicker cells, but the shape of Jph{V) is very similar for cells 
with thicknesses between 55 and 130 nm. This is surprising, 
since e.g. polaron pair dissociation and thereby 7ph are sup- 
posed to be dependent on electric field E. 21 In the simplest ap- 
proximation, assuming a constant slope in energy bands over 
the whole extension of the device, the effective field is the 
fraction of effective voltage and active layer thickness: 



E = \V-V ms \/d 



(2) 



Characteristic points of the photocurrent, indicated by red 
dashed lines in Figure [5] lie on top of each other when plot- 
ted against effective voltage \V — Vpos | ■ When plotted against 
the calculated field E (Equation Q), the curves of 7ph disperse 
(Figure |5|b)). This indicates that the approximation for the in- 
ternal field in Equation (|2]) is not realistic. The reason for this 
is the voltage drop at the contacts (Figure pi, which greatly 
reduces the electric field in the bulk of the cell. In addition the 
photocurrent is not only governed by field dependent polaron 
pair dissociation, but also depends on the extraction of charge 
carriers. As described below, the active layer thickness does 
not influence this extraction mechanism. 



D. Polaron pair dissociation and charge extraction 

To describe the voltage dependence of the photocurrent, we 
propose a model that takes field-dependent polaron pair (PP) 
dissociation and voltage dependent polaron extraction by drift 
and diffusion into account. A well-known model to describe 
field-dependent PP dissociation was presented by Braun in 
1984. 6 Using the field dependence for ion-pair dissociation by 
Onsager and Langevin recombination for polarons, 7 15 Braun 
derives the following PP dissociation probability: 
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FIG. 5. (Color online) Plot of the photocurrent in reverse direction 
(with the origin shifted to POS) against effective voltage (a) and 
calculated field (b) for cells with different active layer thicknesses. 
While the characteristic points, indicated by the red dashed lines, lie 
on top of each other in the left graph, they disperse in the right graph. 



where kd(E) = Kd(E)fi is the PP dissociation rate. It can 
be calculated by detailed balance from the recombination rate 
and the Coulomb binding energy, which depends on the PP ra- 
dius rpp. kf is the decay rate of a PP to the ground state. Since 
/j and Tf = k^ 1 do not influence PBraun(£) independently, we 
use their product /Jif as a single parameter.^ 

Upon successful separation, the polarons need to travel to 
the electrodes to be extracted. This part can be described with 
a term introduced by Sokel and Hughes.^ They solved the 
problem of photoconductivity in insulators analytically, ne- 
glecting dark current, trapping and recombination. The neg- 
ligence of recombination is a resonable approximationPSHSl 
The result depends on temperature T and applied voltage V: 



Jph — Jph,i 



exp(qV/kT + 1) 2kT 
exp(qV/kT-i) ~ IjV 



(4) 



Using an effective voltage V e ff = \V — Vpos I an d a calculated 
field E = \V e ff\/d, and adding the constant offset ./ph(Vpos)> 
the overall photocurrent can be written as: 
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FIG. 6. (Color online) Comparison of experimental photocurrent 
(cell with Ca/Al cathode under 1 sun) and a model based on Braun- 
Onsager theory and extraction. Field dependent part of Jp^ in for- 
ward and reverse direction (left axis) and PP dissociation probabil- 
ity (right axis) vs. effective voltage. The black lines show the PP 
dissociation after Braun (dotted), extraction after Sokel & Hughes 
(dashed) and their combined product (solid). 
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(5) 

This combination of Braun-model for polaron pair disso- 
ciation and the Sokel-Hughes term for charge extraction was 
first proposed by Mihailetchi et al., 9 but with the assumption 
that V e ff = |V — Vb| ■ This corresponds to the special case of 
a vanishing offset, where Vpos equals Vq. Also, they used 
a Gaussian distribution of polaron pair radii rpp to calculate 
the dissociation efficiency in the Braun model. While this as- 
sumption seems reasonable, a fixed value for rpp could de- 
scribe our data better than a distributed one. 

As shown in Figure [6] the product of PP dissociation af- 
ter Braun and extraction after Sokel-Hughes (Equation (|5]l) 



FIG. 7. (Color online) Effective photocurrent /pj, — iph(Vpos) vs - 
applied voltage for two solar cells with Ca/Al and Ag cathode (sym- 
bols). For both cells, the experimental data (symbols) agrees well 
with the model, Equation |5} (solid lines). 



describes the experimental photocurrent well in both forward 
and reverse directions. A set of parameters in a narrow range 
could be employed to describe the experimental photocurrent 
of solar cells with different active layer thickness and cathode 
material over the whole measured range (Figure [7J. 

Using Equation (|2|i to calculate E, effective values of about 
400 nm for d were used in the model, almost independent of 
measured device thickness and clearly exceeding it. This is a 
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consequence of the voltage drops at the contacts (Figure |2j. 
Low injection barriers result in a high band bending in vicin- 
ity of the electrodes and reduced electric field in the bulk of 
the cell, corresponding to a seemingly greater d. The approx- 
imation E = V e ft/d is therefore flawed. 

The saturated, voltage-dependent part of the photocurrent 
■^Ph.max is typically 5-8 mA cm~ 2 , with higher values for 
thicker cells. The dielectric constant £,- was set to 3.5 and 
temperature T to 310 K (slightly above room temperature to 
account for heating under illumination). The f/Zf product of 
spatially averaged mobility /j and PP lifetime (with respect to 
the ground state) T/ = 1 /kf was set to l-3T0 _14 m 2 V -1 and 
the PP radius rpp to 1.8-2.0 nm. 

Using this parameter set, the dissociation efficiency of po- 
laron pairs in the bulk can be calculated to be 40-60 % at 
zero field. However, extraction is the limiting factor at low 
effective voltages. Only at higher effective voltages, the field 
dependence of polaron pair dissociation determines the shape 
of the photocurrent. 

IV. CONCLUSIONS 

By use of a pulsed measurement technique, we find a 
point of symmetry for the photocurrent at 0.52-0.64 V 
in P3HTPCBM organic solar cells, which agrees with re- 
sults from capacitance-voltage measurements. Based on our 
macroscopic simulation, we identify this voltage as the point 
of flat bands in the bulk of the solar cell, the quasi flat band 
case, well below the built-in potential. We find a voltage- 



independent offset of the photocurrent, which is crucial for 
high /sc, which depends both on thermal annealing conditions 
and cathode material. We also propose a possible origin of 
this offset with the band-bending close to the contacts. Active 
layer thickness has an influence on this offset, but surprisingly 
little impact on the shape of 7ph. 

A model including polaron pair dissociation, based on 
Braun-Onsager theory, and extraction, based on the work by 
Sokel & Hughes, can describe our data for different solar cells 
with a narrow range of parameters. The electric field in the 
bulk of the cell is greatly reduced by the injection barriers. 
At zero electric field, polaron pair dissociation efficiency is 
40-60%. 

Since the maximum power point is close to Vpos, where 
the effective voltage is zero, the photocurrent under operating 
conditions is largely determined by the voltage independent 
offset. For best device performance it is necessary to maxi- 
mize this constant contribution from the contact regions and 
simultaneously enhance polaron pair dissociation and polaron 
extraction in the bulk. 
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